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The influence of losses in the interferometric generation 
and the transmission of continuous- variable entangled light is 
studied, with special emphasis on Gaussian states. Based on 
the theory of quantum-state transformation at absorbing di- 
electric devices, the amount of entanglement is quantified by 
means of the relative-entropy measure. Upper bounds of en- 
tanglement and the distance to the set of separable Gaussian 
states are calculated. Compared with the distance measure, 
the bounds can substantially overestimate the entanglement. 
In particular, they do not show the drastic decrease of en- 
tanglement with increasing mean photon number, as does the 
distance measure. 
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I. INTRODUCTION 

Entangled quantum states containing more than one 
photon on average have been of increasing interest (see 
for example [0) for several reasons. So, maximally en- 
tangled continuous-variable states of EPR type require 
infinitely large mean numbers of photons. Another rea- 
son is that such states might be more robust against de- 
coherence. If, for example, in an experiment with Bell 
states one photon is lost, e.g. absorbed, during the trans- 
mission from a sender to a receiver, the entanglement is 
immediately gone (the state is projected onto a separa- 
ble state). States with more than one photon on average 
have the advantage that even a few photons might get 
lost on the way (leaving behind some mixed state), but 
inseparability might still be preserved, that is, maximally 
entangled photon pairs could still be extracted by means 
of purification procedures |2|-|^] . 

The aim of the present article is to investigate entangle- 
ment properties of bipartite quantum states that generi- 
cally live in infinite-dimensional Hilbert spaces. Typical 
examples are Gaussian states such as two-mode squeezed 
vacuum states, which are the states commonly used in 
quantum communication of continuous-variable systems 
. They are also the states we will be looking at in what 
follows. 

In real experiments, both in quantum-state genera- 
tion and processing, e.g., transmission through generi- 
cally lossy optical systems such as fibers, the necessarily 
existing interaction of the light fields with dissipative en- 
vironments spoils the quantum-state purity, leaving be- 
hind statistical mixtures. Unfortunately, quantification 
of entanglement for mixed states in infinite-dimensional 
Hilbert spaces is yet impossible in practice. It typically 



involves minimizations over infinitely many parameters, 
as it is the case for the entropy of formation as well as 
for the distance to the set of all separable quantum states 
measured either by the relative entropy or Bures' metric 
Q. It is, however, possible to derive upper bounds on the 
entanglement content |^] by using the convexity property 
of the relative entropy. For Gaussian states, however, we 
derive an upper bound based on the distance to the set 
of separable Gaussian states which is far better than the 
convexity bound. These bounds may be very useful for 
estimation of the entanglement degradation in real quan- 
tum information systems. 

The article is organized as follows. In Section || the in- 
fluence of losses in the interferometric entanglement gen- 
eration at a beam splitter is studied. A typical situation 
in quantum communication is considered in Section III, 
in which the entanglement degradation of a two-mode 
squeezed vacuum (TMSV) state is transmitted through 
a noisy communication channel say, two lossy optical 
fibers, is examined. Some concluding remarks are given 
in Section IV 



II. ENTANGLEMENT GENERATION BY 
MIXING SQUEEZED VACUA AT A BEAM 
SPLITTER 

A. Lossless beam splitters 

Let us first consider the case of a lossless beam splitter 
and (quasi-)monochromatic light of (mid-)frequency uj 
(Fig. 0). It is well known that a lossless beam 

splitter transforms the operators of the incoming modes 
ai(w) and 0.2 (cj) to the operators of the outgoing modes 
61(0;) and 62 (ti^) according to 



61 (w) 

1)2 (uj) 



ai{uj) 



(1) 



where T(w) is the unitary characteristic transformation 
matrix of the beam splitter. Equivalently, the operators 
can be left unchanged and instead the density operator 
is transformed with the inverse matrix T^^{uj) — T'^{uj) 
according to 



Pout — Q'n 



Tn^)i°'<'"'),T-(„)p;;"! 

\a2{uj)J \al{uj) 



(2) 



Let each of the two incoming modes be prepared in a 
squeezed vacuum state, i.e., 
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^in = |*i„)(*i„| (3) 

where 

|*in) =5i52|0,0), (4) 

with Si (j = l, 2) being the (single-mode) squeeze operator 

5, = exp[-i(e,ar-Ca,2)" 

= exp(-i,,ar) (l-k.n''"'+^^/%xp(i,?a?) (5) 

{qi ~ tanh |^i|e"^', (/>i = argi^^). Here, the second Hne 
foUows from general disentangling theorems ||l^,|l5|. By 
Eq. (H), the output quantum state is 



where 



Pout = |'I'out)(*out|, 



\^out)=[{l-\qi\'){l-\q2\')f' 
X exp -igi(Tiia| +T2i4)^ 

-^52(^124 +7^224)1 |0,0) 



(6) 



(7) 



The Tij are the elements of the characteristic transfor- 
mation matrix T (at chosen mid- frequency) , which can 
be given, without loss of generality, in the form of 



T R 

-R* T* 



(8) 



with T = |T|e'"^^ and R = |i?|e''^« being the (complex) 
transmission and reflection coefficients of the beam split- 
ter. 

From inspection of Eq. (^) it is seen that the prepara- 
tion of an entangled state is controlled by the parameter 



When ^12 = is valid, then the output state is separa- 
ble. This is the case for — (p2 +2{(I)t — 4>r) = and 
|5'i| = |q2|. On the other hand, if again |gi| = |g2| = l^l 
but 01 - (j)2 +2{(j)T -(I>r)=±tt, then for |ri?| = 1/2 the 
output quantum state is just a TMSV state. 



l^-out) = ITMSV) = yr^kF exp 
where 



-qa\al 



|0,0), (10) 



(11) 



Since, according to Eq. (|6|), the output quantum state 
is a pure state, entanglement is uniquely measured by the 
von Neumann entropy of the (reduced) quantum state of 
either of the output modes. 



£^(^out) = 5*1(2) = -Tr[^i(2) In ^1(2)], 



(12) 



where ^1(2) denotes the (reduced) output density opera- 
tor of mode 1(2), which is obtained by tracing ^out with 
respect to mode 2(1). The result of the numerical calcu- 
lation is illustrated in Figs. I and I for a 50%/50% beam 
splitter. In Fig. || the phases are chosen such that the 
relation — 02 +2(0t —4'b.) — is valid, thus leading 
to a TMSV state for |(7i| = |g2|- Figure ^ shows the case 
where 0i — 02 +2(0t — 0_r) = 0, so that for = \q2\ no 
entanglement is observed. Note that using squeezed co- 
herent states instead of squeezed vacuum states does not 
change the entanglement. This is due to the fact that co- 
herent shifts are unitary operations on subsystems which 
leave any entanglement measure invariant. 



B. Lossy beam splitters 

In practice there are always some losses and things get 
slightly more complicated. The SU(2) group transforma- 
tion in Eq. (|l|) has to be replaced by a SU(4) group trans- 
formation, where the unitary transformation acts in the 
product Hilbert space of the field modes and the device 
modes p6|-p^. As a result, Eqs. (|l|) and (||), respectively, 
have to be replaced by 



and 



,(F) 



/3(w) A(w)a(w) (13) 



TY^''^gjA+{u;)a{u;),A^{iu)a\u;)] , (14) 



where the "four- vector" notation a{Lu) for abbreviating 
the list of operators ai(a;), d2{(^), 5i(ijj), 52(w) [and f3{uj) 
accordingly] has been used. The SU(4) group element 
A(cj) is expressed in terms of the characteristic transfor- 
mation and absorption matrices T(ti;) and A.{lu) of the 
beam splitter as 



62 - giTiiTai + g2ri2T22 - -qiTR* + q2RT* . (9) A(w) 



T(w) A{uj) 
'S{u)C-^{uo)T{uj) C{uo)S-^{uj)A{uj) 



(15) 



with the commuting positive Hermitian matrices 



c(c^) = VtRt+R , s(c^) = v/aRa+M- (16) 

From the above, the output density matrix in the Fock 
basis can be given in the form of (Appendix 

(mi,m2|eiutl"i:"-2) 



^ (l~|glP)(l-k2P) 

mi\m2\ni\n2\ 
1 



(-1) 



mi+m2+ni+ri2 



M 



(0)H:^«.....3(0)> (17) 



n5i!52! ' '''''' 

where IlJJ^(O) denotes the Hermite polynomial of four 
variables with zero argument, generated by the symmet- 
ric matrix M with elements 



2 



giAiiAji + g2Ai2A 



J2 



(18) 



Note that in Eq. (p7| ) it is assumed that the device is 
prepared in the ground state. 

In order to quantify the entanglement content of a 



mixed state ^, such as g^^^l in Eq. (17), we make use of 
the relative entropy measuring the distance of the state 
to the set S of all separable states a j|] , 



E{g) — minTr [^( In g — In tr)] 



(19) 



For pure states this measure reduces to the von Neu- 
mann entropy ( p^ ) of either of the subsystems which can 
be computed by means of Schmidt decomposition of the 
continuous- variable state Q . It is also known that when 
the quantum state has the Schmidt form, 



g = Cn,7n I (f>n , V'n) (0m , V'm | , 



(20) 



then the amount of entanglement measured by the rela- 
tive entropy is given by |19,EQl 



E{g) = -^C„^„lnC„,„ - S{g). 



(21) 



Unfortunately, there is no closed solution of Eq. ([l9| ) for 
arbitrary mixed states. Nevertheless, upper bounds on 
the entanglement can be calculated 0, representing the 
quantum state under study in terms of states in Schmidt 
decomposition and using the convexity of the relative en- 
tropy. 



-^(X!^'"^") - ^PnEig,i), ^Pn = 1. 



(22) 



Applying the method to the output quantum state in 
Eq. (|l7|), i.e., rewriting it in the form of 



-(F) 
Sout 



J2 Ckj,o\k,k){l,l\ 

k,l=Q 

oc oo 

+ ^ ^ Ck,Lrn\k + m,k){l + mj\ 

m=l k,l=0 
oo oc 

+ ^ ^ Ck,Lm\k,k + m){l,l + m\ 



m=l k,l=0 

^ PoQo 



ni — 1 

the inequality (|2^) leads to 

oc 

,(F) 



m—1 



(23) 



E{gZl) < PoEiga) + ^ p„, [E{g„,,,) + Eig^^2)] , (24) 

m — l 

where E{go), E{g„i i), and E(gm.2) can be determined 
according to Eq. (|2^). In the numerical calculation we 



have used the dielectric-plate model of a beam splitter, 
taking the T and A matrices from |l^,^ . The result is 
illustrated in Fig. |^, which shows the dependence on the 
plate thickness of the upper bound of the attainable en- 
tanglement. The oscillations are due to phase matching 
and phase mismatch at certain beam splitter thicknesses 
[cf . Eq. (^)] . Note that the local minima of the curve for 
the lossy beam splitter never go down to zero as do the 
corresponding minima of the curve for the lossless beam 
splitter. This obviously reflects the fact that the result 
for the lossless beam splitter is exact, whereas that for 
the lossy beam splitter is only an upper bound. 



III. ENTANGLEMENT DEGRADATION IN 
TMSV TRANSMISSION THROUGH LOSSY 
OPTICAL FIBERS 

Let us now turn to the problem of entanglement degra- 
dation in transmission of light prepared in a TMSV 
state through absorbing fibers. The situation is some- 
what different from that in the previous section, since 
we are effectively dealing with an eight-port device as 
depicted in Fig. H, where the two channels are character- 
ized by the transmission (T^) and reflection [Ri) coeffi- 
cients (j = l,2). In particular for perfect input coupling 
(i?i = 0), the system is essentially characterized by the 
transmission coefhcients T^. 

From Eq. (|l^) it is easily seen that in the Fock basis a 
TMSV state reads 

OO 

ITMSV) = yT~kf E(-'?)"l"' (25) 

n=0 

whose entanglement content is 

£;(|TMSV)) = - In {l-\q\^) In \q\\ (26) 

Application of the quantum-state transformation (|lj) 
yields (i?, = 0) 



OO OO 



-(F) 
fi'out 



m=0 fe,i=0 

X (c™|m+A:)(/c| +H.C.) (g) ((i™|m-|-/)(;| -l-H.c.) , (27) 



where 



and 



k.l.rn 



-g)"/2r™ (1 - iJ™o) , 

-Q)"/2r™ (1 - ij^o) , 

[|gp(l-|Tip)(l-|r2n]V(a+m)! 



(28) 
(29) 



^fc!Z!(fc+TO)!(? + m)!(a-fc)!(a-/)! 



1^1 1 



l-lTiP 



IT, 
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(30) 



[a = max(fc, Z)]. Note that in Eq. ( pTj) the fibers are as- 
sumed to be in the ground state. 



A. Entanglement estimate by pure state extraction 

The amount of entanglement contained in the (mixed) 
output state (p^ can also be estimated, following the 



line sketched in Section II B. In particular, the convexity 



of the relative entropy can be combined with Schmidt 
decompositions of the output state in order to calculate, 
on using the theorem (pl|), appropriate bounds on entan- 
glement. Before doing so, let us first consider the case of 
low initial squeezing, for which the entanglement can be 
estimated rather simply. 



Note that for Ti =T2 = 1 the entanglement of the TMSV 
state is preserved, i.e., Eq. (^) reduces to Eq. (^. In 
Fig. ^, the estimate of entanglement as given by Eq. ( |3^ ) 
is plotted as a function of the transmission length and 
the strength of initial squeezing for Ti ~T2—T, where T 
is given by the Lambert-Beer law of extinction, 



l/lj, 



(36) 



Here, is the real part of the complex refractive in- 
dex, Ia —c/{niuj) is the absorption length, and I is the 
transmission length. 

It is worth repeating that the estimate given by 
Eq. (33) is valid for low squeezing only. Higher squeez- 
ing amounts to more excited density matrix elements and 
Eq. (|3^) might become wrong. Moreover, we cannot even 
infer it to be a bound in any sense since no inequality has 
been involved. A possible way out would be to extract 
successively more and more pure states from Eq. (|2^). 
But instead, let us turn to the Schmidt decomposition. 



1. Extraction of a single pure state 



2. Upper bound of entanglement 



Since, by Eqs. ( pT] ) - (|3C|), for low squeezing only a few 
matrix elements are excited which were not contained in 



the original Fock expansion (25), we can forget about the 
entanglement that could be present in the newly excited 
elements and treat them as contributions to the separable 
states only. Following [0 , the inseparable state relevant 
for entanglement can then be estimated to be the pure 
state 



yi~A|*) 



'1 



000 



KoonCndn\n, n) . (31) 



It has the properties that only matrix elements of the 
same type as in the input TMSV state occur and the 
coefiicients of the matrix elements |0, 0) |n, n) are met 
exactly, i.e., 

(l-A)(0,0|vI')(vl/|n,n) - {G^QolW^n). (32) 

In this approximation, the calculation of the entangle- 
ment of the mixed output quantum state reduces to the 
determination of the entanglement of a pure state |17[: 



{i-x)Em) 



l-x 



In 



-y 



l-x 



l-x) 



-y 



where 



{l-x){[y+{l-x)^] Hl-x)-ylny} 
[y-(l-x)2]2 

x=\qf{l-\T,f){l-\T2f), 



y 



(33) 



(34) 



(35) 



In a similar way as in Section [IB, an upper bound 
on the entanglement can be obtained 0, if the density 
operator (^^ is rewritten as the convex sum of density 
operators in Schmidt decomposition. 



£ul^ E \Ak,i\k,k){l,l\ 

k,l^O ^ 
oo 

+ ^ Bkj,m\k-\-m,k){l-\-mJ\ 

m— 1 

oo V 

+ ^Ck,i,m\k,k + m){l,l + m\y (37) 

m=l ^ 



and the inequality (|2^) together with Eq. ( 21 ) is applied. 
The result is illustrated in Fig. 0. 

From general arguments one would expect the entan- 
glement to decrease faster the more squeezing one puts 
into the TMSV, because stronger squeezing is equivalent 
to saying the state is more macroscopically non-classical 
and quantum correlations should be destroyed faster. As 
an example, one would have to look at the entangle- 
ment degradation of an n- photon Bell- type state l^*^), 
£^(|*^)) < |T|2"ln2 ||. Since the transmission coef- 
ficient T decreases exponentially with the transmission 
length, entanglement decreases even faster. Note that 
similar arguments also hold for the destruction of the in- 
terference pattern of a cat-like state ^ ja) -|- | — a) when 
it is transmitted, e.g., through a beam splitter. It is well 
known that the two peaks (in the jth output channel) de- 
cay as \Tji p, whereas the quantum interference decays as 



|r,-ipexp[-2|ap(l- 



The upper bound on the entanglement as calculated 
above seems to suggest that the entanglement degrada- 
tion is simply exponential with the transmission length 
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for essentially all (initial) squeezing parameters, which 
would make the TMSV a good candidate for a robust en- 
tangled quantum state. But this is a fallacy. The higher 
the squeezing, the more density matrix elements are ex- 
cited, and the more terms appear, according to Eq. (37), 
in the convex sum (|2^). Equivalently, more and more 
separable states are mixed into the full quantum state. 
By that, the inequality gets more inadequate. In order to 
see this better, we have shown in Fig. ^ the upper bound 
on the entanglement for just two different (initial) squeez- 
ing parameters \q\ =0.71 (equivalent to the mean photon 
number of fi = 1, solid line) and |g| = 0.9535 (n = 10, 
dashed line). For small transmission lengths, hence very 
few separable states are mixed in, the curves show the 
expected behavior in the sense that the state with higher 
initial squeezing decoheres fastest. The behavior changes 
for larger transmission lengths. We would thus conclude 
that the upper bound proposed in is insufficient. 



B. Distance to separable Gaussian states 

The methods of computing entanglement estimates 
and bounds as considered in the preceding sections are 
based on Fock-state expansions. In practice they are typ- 
ically restricted to situations where only a few quanta of 
the overall system (consisting of the field and the device) 
are excited, otherwise the calculation even of the matrix 
elements becomes arduous. Here we will focus on an- 
other way of computing the relative entropy, which will 
also enable us to give an essentially better bound on the 
entanglement (for other quantities that characterize, in 
a sense, entanglement, see [p3[). 

Since it is close to impossible to compute the distance 
of a Gaussian state to the set of all separable states we re- 
strict ourselves to separable Gaussian states. A quantum 
state is commonly called Gaussian if its quantum char- 
acteristic function is Gaussian. By the general relation 
for a iV-mode quantum state 



1 



d^'^ax(-a)£'(a) 



(38) 



it is obvious that the density operator of a Gaussian state 
can be written in exponential form of 



■ A/'exp 



a) Ma 



(39) 



where Mo- is a Hermitian matrix that can be assumed 
to give a symmetrically ordered density operator, and J\f 
is a suitable normalization factor. Here and in the fol- 
lowing we restrict ourselves to Gaussian states with zero 
mean. Since coherent displacements, being local unitary 
transformations, do not influence entanglement, they can 
be disregarded. 

The relative entropy (M) can now be written as 



£;^(g) = minTr |p In p - InW + (a^ a) ^ " ^ | 

= Tr (^ In ^) + rnin {a^ a) ^ - In . (40) 

Since we have chosen the density operator a to be sym- 
metrically ordered, the last term in Eq. (^0|) is nothing 
but a sum of (weighted) symmetrically ordered expecta- 
tion values {a^^a"-) s=o {m + n = 2). For a Gaussian quan- 
tum state g it can be shown (Appendix ^) that Eq. ( |40| ) 
can equivalently be written in terms of the matrix Dg in 
the exponential of the characteristic function of g as 



ER{g) = Tr(^ln^) 



min[iTr (M^D^) - InW] . (41) 



From the above it is clear that we only need the matrix 
T)g (which is unitarily equivalent to the variance matrix). 
For a Gaussian distribution with zero mean the elements 
of the variance matrix V are defined by Vij = {(i, Cj)s=o 
as the (symmetrically ordered) expectation values of the 
quadrature components C = (ii,Pi: X2,P2)- 

The variance matrix of the TMSV state (25) reads 
(g = tanh|e|e^^e=|e|e^*) 



X Z 

ZT Y 



(42) 

with the notation c = cosh2|^|, si = sinh 2|^| cos (/>, and 
S2 — sinh 2|^| sin <f>. In case = the variance matrix (|42 
reduces to the generic form 




Vn = 



X zi \ 

X Z2 

zi y 

V Z2 y 



(43) 



specified by four real parameters. Note that the variance 
matrix of any Gaussian state can be brought to the form 
(H) by local Sp(2,R)(g)Sp(2,M) transformations 0, so 
that we can restrict further discussions to that case. 

Application of the input-output relations ( p^ gives for 
the elements of the variance matrix of the output state, 
on assuming that the two modes are transmitted through 
two four-port devices prepared in thermal states of mean 
photon numbers rithi, [^ | 

All = A22 = 

iclTip + i|i?i|2 + (nthi + i) (l-|ri|2-|i?ip) , (44) 
Yii = Y22 = 

\C\T2? + i|i?2p + K2 + i) {l-\T2? - \R2?) , (45) 



Zii = -Z22 = -\sRc{TiT2), 
Z\2 = Z2\ = -5sIm(TiT2) 



(46) 
(47) 
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((/> = 0). With regard to optical fibers with perfect in- 
put coupUng {Ri = 0) and equal transmission lengths, we 
again may set \Ti\ = e^'^^^ . Moreover, we may assume 
real Ti and thus set Z12 — Z21 ~ 0. 

First, one can check for separability according to the 
criterion [||,|4| 

detXdetY+ (i - |detZ|)^ - Tr (^XJZJYJZ'^j) 

> i (det X + det Y) , (48) 
which reduces to 



A{xy - zt){xy - z^) > {x^ + y^) + 2|ziZ2| 



(49) 



Combining Eqs. (jij) - (^Qj), it is not difficult to prove 
that the boundary between separability and inseparabil- 
ity is reached for P, |l7| , p5t 



I = h 



Ia 



In 



1 



-2|?| 



(50) 



It is worth noting that this is exactly the same con- 
dition as for the transmitted state still being a squeezed 
state or not. To show this, we calculate the normally- 
ordered variance (: (A_F)^) :) of a phase-sensitive quan- 
tity such as F= iFile^'^ifli + |F2|e'"^=^a2 -I- H.c. Using the 
input-output relations (p^, the normally ordered vari- 
ance of the output field is derived to be 

(: (AF)2) = 2\Fif [|Ti f sinh^ \C\+nthi (l - \Tif)] 
+ mf [|T2|'sinh2 +nth2 (1 - \T2f)] 
-2\FiF2TiT2\smh2\^\ cos{ipi+ip2 + ipT + 4>) (51) 

[Tj = I I e"^^' , i — l,2; (pr — VTi + </5T2 ] ■ For equal ampli- 
tudes |Fi| = IF2I = |F| and equal fibers |Ti| = IT2I = |T|, 
Hthi = nth2 = nth the (phase-dependent) minimum is ob- 
tained to be 



:(AF)2):)out 



4|^^|' 



nth(l-|T|2)-|r|2sinh|e|£ 



(52) 



Equation (B3) exactly leads to the condition (pQ), i.e 



(:(AF)2):)„,, 



< if l<ls 
> if l>ls 



(53) 



Therefore, measurement of squeezing corresponds, in 
some sense, to an entanglement measurement. 

In order to obtain (for l<ls) a measure of the entangle- 
ment degradation, we compute the distance of the output 
quantum state to the set of all Gaussian states satisfying 
the equality in (|4^), since they just represent the bound- 
ary between separability and inseparability. These states 
are completely specified by only three real parameters 
[one of the parameters in the equality in (^) can be com- 
puted by the other three]. With regard to Eq. (41), min- 
imization is thus only performed in a three-dimensional 



parameter space. Results of our numerical analysis are 
shown in Fig. ^. It is clearly seen that the entangle- 
ment content (relative to the entanglement in the initial 
TMSV) decreases noticeably faster for larger squeezing, 
or equivalently, for higher mean photon number [the rela- 
tion between the mean photon number n and the squeez- 
ing parameters being n = sinh^ |^| = |(7|^/(1 — 

It is very instructive to know how much entanglement 
is available after transmission of the TMSV through the 
fibers. Examples of the (maximally) available entan- 
glement for different transmission lengths are shown in 
Fig. 10, One observes that a chosen transmission length 



allows only for transport of a certain amount of entangle- 
ment. The saturation value, which is quite independent 
of the value of the input entanglement, drastically de- 
crease with increasing transmission length (compare the 
upper curve with the two lower curves in the figure). This 
has dramatic consequences for applications in quantum 
information processing such as continuous-variable tclc- 
portation, where a highly squeezed TMSV is required in 
order to teleport an arbitrary quantum state with suffi- 
ciently high fidelity ^ . Even if the input TMSV would be 
infinitely squeezed, the available (low) saturation value of 
entanglement principally prevents one from high-fidelity 
teleportation of arbitrary quantum states over finite dis- 
tances. 



C. Comparison of the methods 

In Fig. prl th e entanglement degradation as calculated 
in Sectionlni B is compared with the estimate obtained in 



Section III A 1 and the bound obtained in Section III A 2 



The figure reveals that the distance of the output state 
to the separable Gaussian states (lower curve) is much 
smaller than it might be expected from the bound on 
the entanglement (upper curve) calculated according to 
Eq. (H) together with Eqs. (|lj) and (|^), as weU as the 
estimate (middle curve) derived by extracting a single 
pure state according to Eq. (BSh. Note that the entan- 
glement of the single pure state (|3^) comes closest to the 
distance of the actual state to the separable Gaussian 
states, whereas the convex sum ( |37| ) of density operators 
in Schmidt decomposition can give much higher values. 
Both methods, however, overestimate the entanglement. 
Since with increasing mean photon number the convex 
sum contains more and more terms, the bound gets worse 
[and substantially slower on the computer, whereas com- 
putation of the distance measure ( pl| ) does not depend 
on it]. 

Thus, in our view, the distance to the separable Gaus- 
sian states should be the measure of choice for determin- 
ing the entanglement degradation of entangled Gaussian 
states. Nevertheless, it should be pointed out that the 
distance to separable Gaussian states has been consid- 
ered and not the distance to all separable states. We 
have no proof yet, that there does not exist an insepa- 
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rable non-Gaussian state which is closer than the closest 
Gaussian state. 



IV. CONCLUSIONS 

In the present article the interferometric generation 
and the transmission of entangled light have been stud- 
ied, with special emphasis on Gaussian states. The opti- 
cal devices such as beam splitters and fibers are regarded 
as being dispersing and absorbing dielectric four-port de- 
vices as typically used in practice. In particular, their 
action on light is described in terms of the experimen- 
tally measurable transmission, reflection, and absorption 
coefficients. 

An entangled two-mode state can be generated by mix- 
ing single-mode non-classical light at a beam splitter. 
Depending on the phases of the impinging light beams 
and the beam-splitter transformation, the amount of en- 
tanglement contained in the outgoing light can be con- 
trolled. For squeezed vacuum input states and appropri- 
ately chosen phases, maximal entanglement is obtained 
for lossless, symmetrical beam splitters. In realistic ex- 
periments, however, losses such as material absorption 
prevents one from realizing that value. 

When entangled light is transmitted through optical 
devices, losses give always rise to entanglement degrada- 
tion. In particular, after propagation of the two modes of 
a two-mode squeezed vacuum through fibers the available 
entanglement can be drastically reduced. Unfortunately, 
quantifying entanglement of mixed states in an infinite- 
dimensional Hilbert spaces has been close to impossible. 
Therefore, estimates and upper bounds for the entangle- 
ment content have been developed. 

The analytical estimate employed in this article is 
based on extraction of a single pure state from the out- 
put Gaussian state, using its reduced von Neumann en- 
tropy as an estimate for the entanglement. However, this 
method is neither unique, since there are many different 
ways of extracting pure states, nor is it an upper bound, 
since nothing is said about the residual entanglement 
contained in the state which is left over. In principle, 
more and more pure states could be extracted until the 
residual state becomes separable. 

Instead, an upper bound can be calculated by decom- 
posing the output Gaussian state in a convex sum of 
Schmidt states as proposed in 0. The disadvantage of 
this method is that the bound gets worse for increasing 
(statistical) mixing. In particular, it may give hints for 
large entanglement even if the quantum state under con- 
sideration is almost separable. 

In order to overcome the disadvantage, the distance of 
the output Gaussian state to the set of separable Gaus- 
sian states measured by the relative entropy is consid- 
ered. It has the advantage that separable states obvi- 
ously correspond to zero distance. Although one has yet 
no proof that there does not exist a non-Gaussian sep- 



arable state which is closer to the Gaussian state under 
consideration than the closest separable Gaussian state, 
one has good reason to think that it is even an entangle- 
ment measure. In any case, it is a much better bound 
than the one obtained by convexity. In particular, it 
clearly demonstrates the drastic decrease of entanglement 
of the output state with increasing entanglement of the 
input state. Moreover, one observes saturation of entan- 
glement transfer; that is, the amount of entanglement 
that can maximally be contained in the output state is 
solely determined by the transmission length and does 
not depend on the amount of entanglement contained in 
the input state. 
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where the Hermite polynomials of four variables are gen- 
erated by the symmetric matrix M with elements 



M,. 



(A5) 



Using the relation between Hermite polynomials of one 
variable and those of several variables |Eq| , 
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where 



APPENDIX A: FOCK-STATE EXPANSION OF 
MULTIMODE SQUEEZED VACUUM STATES 

Let us consider an incoming field prepared in the 
squeezed vacuum state (^ and an absorbing beam split- 
ter in the ground state. The quantum-state transforma- 
tion formula (|H) then leads to 



fout 



E (.9i,52|^a;(6)^ai(e2)|0,0,0,0) 

31,92=0 

(0,0,0,0|5^,(6)^L(a)l5i,52) (Al) 



where the transformed operators a[ are defined by 



EA, 



(A2) 



according to the rules of quantum-state transformation. 
Equivalently, the Fock-state expansion of the density ma- 
trix reads 

{mi,m2\gZ't\ni,n2) = 
(mi, m2|(5i,ff2|5„-,(6)5„^ (6)10,0,0,0) 
(0,0,0,0|5^, {^2)Sl, (ei)|5i,52)|ni,n2). 



(A3) 



Expanding the Fock states in terms of coherent states 
and using the squeeze operator in the form given in the 
second line in Eq. (^), we obtain after performing all 
integrals 



(a,x) 



hi 



(A7) 



and '/>(a) = (/)(a, a) , we get for the multivariable Hermite 
polynomial of zero argument 
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2/2 terms 



(to = TOi + TO2 + gi + (?2)- Here, the 7^-sum runs over 
all m!/(TOi!TO2!(7i!52!) possible combinations to distribute 
nil mdices 1, TO2 indices 2, gi indices 3, and 92 mdices 4 
among the indices ii, ji, .. ., i„/2, jm/2- 

In par ticular, if we restrict ourselves to two dimensions, 
Eq. ( |A^ ) simplifies to 



^mi, 7772(0) 



^777(0) 

2777/2 



Mr; 



(7771 — 7n2 — 7^)/4 ^^^(7712 — mi — 1/)/4 
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ItooII — II 



77= it/ 



f n\{^J. - 2ny.{n - ij/)! 



M^j^{2Mi2T'^"M22 (A9) 



[n — max(TOi,TO2), V— |toi — TO2I]. The sum can also be 
calculated leading to Gegenbauer, Jacobi or associated 
Legendre polynomials ||27| . 

Another way of writing is the one we used for the nu- 
merical calculation of the density matrix elements. The 
method, however, is only applicable in cases where the 
number of variables the Hermite polynomial depends on 
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is sufBciently small. The multivariable Hermite polyno- 
mial (now in four variables) of zero argument can be writ- 
ten as 



(-1) 



mi+m2+3i+g2 ijM 



Qmi+m2+gi+g2 



mi,m2, 91,92 



(0) = 



exp 



A=o' 



(AlO) 



with M b eing given by Eq. (A5). Expanding the rhs of 
Eq. (AlO), the only surviving term is the one propor- 
tional to (A'^MA)('"i+'"2+9i+92)/2^ Multinomial expan- 
sion of this term then leads to Eq. (p7|). 
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APPENDIX B: MULTIMODE GAUSSIAN 
DENSITY OPERATORS AND WIGNER 
FUNCTIONS 

Given the Wigner function of an A^-mode Gaussian 
state of the form of 



Wn{C) 



1 



(27r)JvVdetV 



exp(-iC^V-ic) , (Bl) 



where C = {xi,pi, . . . ,X]S! ,Pn) is the 2Af-dimensional 
"vector" of the quadrature components of the N (com- 
plex) variables at and V is the 2N x 2N variance matrix 
of the quadrature components. The characteristic func- 
tion defined by the Fourier transform reads 



XN{ri) = exp(-i?7'^Vr7) 



(B2) 



Alternatively, the quantum state can can be given by the 
density operator 



Q = 



exp 


-i(ata)M( 








exp 


- (at a) M 




]} 



(B3) 



with an obvious definition of the matrix D, thus estab- 
lishing a relation between the matrix M in the exponen- 
tial of the density operator and the matrix D in the ex- 
ponential of the characteristic function. From the third 
to the fourth line in Eq. (B6) we have used the expression 
for the characteristic function of a thermal state |^^. In 
due course, the normalization of the density operator is 
obtained as 



N 



TV = J]^ 2 sinh 



(B7) 



The above description shows a way to compute the 
entropy of a Gaussian quantum state g as well as the 
relative entropy between two Gaussian quantum states g 
and a as 



N 



Tr (e In ^) = ^ In ( 2 sinh ^ ) - ^ Tr (M,D,) , (B8) 



N 



i=l 



TT{g\na) = ^ ln( 2 sinh ^"j - 1 Tr (M^D,) , (B9) 



where the Qi and -di are respectively the eigenvalues of 



Mg and 



(where a is actually an A^-dimensional "vector" with 
"components" a^). 

In order to relate the matrix M to the matrix V, we 
introduce a unitary transformation 



u 



(B4) 



where the matrix U is chosen such that it diagonalizes 
M, hence U+MU = (with being diagonal). Note, 
that U satisfies the generalized unitary relation 



UJU^ 



with J = diag (Ijv, —I 



(B5) 



The n, the characteristic function of the density operator 
(H is 
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FIG. 1. Squeezed states and \^2) impinging on a 
beam splitter producing entangled light beams. 
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FIG. 2. Entanglement produced at a lossless 50%/50% 
beam splitter by mixing two modes prepared in squeezed vac- 
uum states as a function of |gi | and \q2 \ for the phase condition 




FIG. 3. Entanglement produced at a lossless 50%/50% 
beam splitter by mixing two modes prepared in squeezed vac- 
uum states as a function of |gi | and 1 for the phase condition 
2((?!'ii — (?!)t) + 02 —01 = 0. 




FIG. 4. Entanglement produced at a lossless beam splitter 
with refractive index n—lAl (dashed curve) as a function of 
the beam splitter thickness I. The full curve shows the upper 
bound of the produced entanglement at a lossy beam splitter 
with n = 1.41 -f O.li. The squeezing parameters chosen are 
qi =q2 =0.5. 
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FIG. 5. 
tum state 
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A two-mode input field prepared in the quan- 
I is transmitted through two absorbing dielectric 
four-port devices, ai, as {a'2, a'4) being the photonic operators 
of the relevant input (output) modes. 
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FIG. 6. Estimate of the entanglement, Eq. (|33[), ob- 
served after transmission of a TMSV through absorbing fibers 
(Ti =T2) as a function of the squeezing parameter \q\^ and the 
transmission length I. 
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FIG. 7. Upper bound on the entanglement degradation of 
a TMSV transmitted through absorbing fibers (Ti =T2) as a 
function of the squeezing parameter \q\ and the transmission 
length I. In the numerical calculation, at most 30 photons per 
mode have been taken into account which is obviously not 
sufficient for higher squeezing when higher photon-number 
states are excited. 
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FIG. 8. Upper bound of the entanglement degradation 
of a TMSV transmitted through absorbing fibers (Ti=T2) 
as a function of the transmission length I for the squeezing 
parameters 5 = 0.1 (solid line) and q — 0.9 (dashed line). In 
the numerical calculation, Fock states \n) up to < 0.02 have 
been taken into account. 
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FIG. 9. Entanglement degradation of a TMSV transmit- 
ted through absorbing fibers {Ti=T2) as a function of the 
transmission length / for the (initial) mean photon numbers 
n= 1 (|g|~0.7071) (topmost curve), fi=10 (|g| ~ 0.9535), 
n= 10^ (|g| ~0.9950), and n=10^ (|g|~0.9995) (lowest 
curve) . 
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FIG. 10. Available entanglement after transmission of a 
TMSV through absorbing fibers {Ti=T2) as a function of 
the squeezing parameter ^ for various transmission lengths / 
[1 = (topmost curve), I = 10"^^^ (middle curve), / = 10~^Ia 
(lowest curve)]. For |C|<0.5 and 1/Ia '^W~^ , the numerical 
accuracy of the values of -B(piul) decreases due to low accu- 
racy in the eigenvector computation. 




FIG. 11. Comparison of the upper bound on entanglement 
(upper curve) according to Fig. [? , the entanglement estimate 
(middle curve) according to Fig. q, and the distance measure 
(lower curve) according to Fig. M for the mean photon number 
n= 1 (|g|~0.7071). 
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